Surprises in the phase diagram of an Anderson impurity model for a single Cgg molecule 
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We find by Wilson numerical renonnalization group and conformal field theory that a three-orbital Anderson 
impurity model for a CgJ~ molecule has a very rich phase diagram which includes non-Fermi-liquid stable and 
unstable fixed points with interesting properties, most notably high sensitivity to doping n. We discuss the 
implications of our results to the conductance behavior of Cso-based single-molecule transistor devices. 
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The characteristic behavior of an Anderson impurity model 
(AIM) emerges often unexpectedly in physical contexts which 
are apparently faraway from magnetic alloys. The typical ex- 
ample is the conductance behavior of nano-scale devices, e.g. 
quantum dots and single-molecule transistors (SMT), where 
Kondo-assisted tunneling may lead to nearly perfect transmis- 
sion at zero bias in spite of a large charging energy 1 1). Usu- 
ally the two conduction leads are bridged by a single quan- 
tum dot or molecular level, which therefore behaves effec- 
tively as a one-orbital AIM. Less common is the case when 
several levels happen to be nearby degenerate thus realizing 
in practice a multi-orbital AIM. This may be achieved for 
instance in SMTs built with high-symmetry molecules. Re- 
cently a zero-bias anomaly has been reported for a Cgo based 
which has been proven to be of the Kondo-type by 
its splitting under the action of a magnetic field or magnetic 
leadsOl- FuUerene's large electron affinity is known to yield 
to electron transfer into Cgo when adsorbed on metallic sub- 
strates. The actual number of doped electrons depends on the 
substrate IJJ, but may also be controlled by attaching alkali 
atoms to the molecule 0]. In this case a variety of conduc- 
tance behavior has been reported depending on the number of 
K-atoms attached to Ceo 1 5], ranging from Kondo-like reso- 
nances to Fano-like anti-resonances. 

Motivated by this opportunity, in this Letter we study the 
behavior of an AIM for Crq" by Wilson Numerical Renor- 
malization Group (NRG)|6] and Conformal Field Theory 
(CFT)0]. We obtain a phase diagram which includes Fermi- 
and non-Fermi-liquid phases, with a doping dependence qual- 
itatively in agreement with experiments. In particular we find 
a Kondo-like zero-bias anomaly for dopiiig n = 1, likely the 
case of Ceo on AufT], as found in Refs.fz, 's']. For n = 2, 
which should correspond to pure or slightly K-doped Ceo on 
Ag0], we predict instead a conductance-minimum at zero- 
bias, compatible with actual observations fs*]. Finally, for 
n = 3 we find a conductance minimum with a non-analytic 
voltage behavior G{V) - G(0) - \V\^/^. 

The LUMO's of Ceo are three-fold degenerate ti„ 
orbitals|8]. The electron-electron interaction acts as if these 
orbitals were effectively p-orbitals. Therefore the Coulomb 
energy of a Cg^," molecule with n valence electrons into a 
configuration with total spin S and angular momentum L is 



Wmt =nu + Hj, with Hu ^ U {n- tiq) /2 and n.j = 
-J [2S{S + 1) + ^L{L + 1)] . The reference valency, no, is 
controlled in SMTs by the gate voltage, the metal substrate 
and the alkali doping. The Coulomb exchange, J > 0, favor- 
ing high-degeneracy states (conventional Hund's rules), com- 
petes in Ceo with the Jahn-Teller coupling to eight fivefold- 
degenerate vibrational modes of Hg symmetry! 8], which, on 
the contrary, prefers low-degeneracy configurations (inverted 
Hund's rules). If the vibrational frequencies, ranging from 35 
to 196 meVfS], are larger than the Kondo temperature, one 
can neglect retardation effects, which amounts to renormalize 
J -> J - SEjt/4:, where Ejt ^ 169 meV is the Jahn-Teller 
energy gain. Although first principle calculations! 9] do pre- 
dict an inverted effective-exchange J — 3Ejt /4 ~ —50 meV, 
in what follows we will consider both J > and J < cases. 

The impurity is coupled to a bath of conduction elec- 
trons, and both the bath and the hybridization are for 
simplicity assumed to be SU(6) and particle-hole (p-h) 
invariant. When discussing the stability of the fixed 
points, we shall take into account deviations from this 
high-symmetry case. The non-interacting Hamiltonian, 
Ho = Hbath + Hhyb can always be written as a one- 
dimensional chain of non-interacting electrons, Hbath ~ 
EamE^>lU[4„^^c^^^„^^ + H.c.y hybridized at one 

edge to the impurity, Hhyb J^am (^m^ + iJ.c.) . 

Here c|,„^ and c - are, respectively, the creation and annihi- 
lation conduction-electron operators at chain site i with spin 
a =1,1 and angular momentum component m = —1,0,1, 
while dl^^ and d„^^ are the analogous impurity operators. 

Although we analyzed by NRG the full AIM, Haim = 
"Hq + Hint, the CFT analysis is more conveniently done for 
the Kondo model onto which the AIM maps for large U. One 
finds that Haim Hbath + H.j + Hk, where the Kondo 
exchange is given by 
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(-l)"'"i/2 (dl^^d^^,^, - 1/2 Smm' Saa'^j is the impurity /7- 
h operator with angular momentum / = 0,1,2 and z- 
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component A, as well as spin s = 0, 1 and z-component 
a, and analogously for the operator at site 1 of the 
chain. C^f'-,-, are Clebsch- Jordan coefficients. Notice 

that {^S ® d ^^^^^ = {no — 3)/V6 is zero when uq = 
3. Since the U(6) symmetry of Hbath is reduced to 
charge U(l), spin SU(2) and orbital 0(3), the proper con- 
formal embedding of the conduction electrons is U(6)i D 
U(l)xSU(2)3xSU(2)8xZ3, where SU(2)3 refers to spin and 
SU(2)8 to the angular momentum.! 10| The 3-state Potts sec- 
tor Za, which reflects orbital permutational-symmetry, can be 
interpreted as emerging out of the charge iso-spin SU(2)3 D 
U(l)xZ3,|7] with generators 4 = - 3)/2, /+ = 

EL-i(-ircLT^l-™i' ^"d/- = (/+)^ 

They commute with Hbath and also with the Kondo-exchange 
terms with {l,s) = (1, 0), (0, 1), (2, 1), while they do not with 
those having (/, s) = (0, 0), (2, 0), (1, 1). We recall that the 
Z3 Potts model has primary fields /, a, Z and e, with dimen- 
sions 0, 1/15, 2/3 and 2/5, respectively, fr'l 

Since the AIM is invariant upon uq — > 6 — rig, we only need 
to consider uq — 1, 2, 3. For no = 1, Hj is ineffective, lead- 
ing at large J7 to a conventional Kondo-screened SU(6) model. 
Conversely, nontrivial behavior may appear for no = 2, 3. 
Here the Kondo screening, controlled by the Kondo tempera- 
ture Tk, is hampered by the exchange splitting J. While the 
former takes advantage from the impurity coherently tunnel- 
ing among all available ng-electron configurations, the latter 
tends to lock the impurity into a sub-set of states with well 
defined S and L. When | J| ~ Tk, neither prevails and a 
nontrivial behavior may emerge. 

In the p-h symmetric case, no = 3, the avail- 
able impurity states have quantum numbers {S,L) — 
(3/2, 0), (1/2, 2), (1/2, 1). For conventional Hund's rules, 
J > 0, the lowest energy state has {S,L) = (3/2,0), while 
for J < {S, L) — (1/2, 1) is favored. Let us first assume 
a positive and very large J ^ Tk- In this case the impu- 
rity effectively behaves as a spin-3/2. If we project Q onto 
the above impurity configuration, only (/, s) = (0, 1) sur- 
vives, which is just the spin-exchange {8V'^/3U) S ■ Si. The 
model reduces to a spin-3/2 impurity coupled to three conduc- 
tion channels, hence it describes a perfectly Kondo-screened 
fixed point which has Fermi-liquid (FL) behavior The im- 
purity density of states (DOS), p{e), is uneffected by T-Cmt at 
the chemical potential, p{Q) — po, meaning a perfect trans- 
mission in the SMT, G/{2e^/h) = p{Q)/po = 1, with G the 
zero-bias conductance per orbital. 

In the opposite extreme of a strong inverted Hund's rule, 
J < -Tk < 0, the impurity locks into the (5, L) = (1/2, 1) 
state. The Kondo exchange projected onto this subspace con- 
tains (/, s) = (1, 0), (0, 1), (2, 1), hence still commutes with 
the iso-spin generators. Thus we expect the asymptotic spec- 
trum to be described by conformal towers identified by the 
quantum numbers of the iso-spin, /, spin, S, and angular mo- 
mentum, L. To infer the fixed-point, let us assume a Kondo 
exchange Jk ^ V'^ /U much larger than the bath bandwidth. 
In this case we have first to solve the two-site problem which 



TABLE I: Left Table: low-energy spectrum at the non-Fermi liquid 
stable fixed point for no = 3, in units of the fundamental level spac- 
ing, as predicted by CFT, Ecft, and as obtained by NRG, Enrg- 
Right Table: dimension x of the most relevant boundary operators. 
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includes the impurity and the site 1 of the chain. The lowest 
energy states are obtained by coupling into an overall singlet 
the impurity with a {S, L) — (1/2, 1) electron-configuration. 
Yet one can form three states at site 1 with {S, L) = (1/2, 1), 
namely through 1, 3 or 5 electrons. The net result is that the 
two-site ground state is threefold degenerate and represents 
an effective iso-spin- 1 impurity. Projecting onto this mani- 
fold the hopping term connecting site 1 to site 2, we get a 
new Kondo exchange acting only among iso-spins, namely 

[til Jk) I ■ l2- This model describes an iso-spin-1 three- 
channel Kondo model which is non-Fermi liquid (NFL). By 
applying the fusion hvpothesis lUlll . which is a formal way in 
CFT to make the impurity dissolve, with its quantum numbers, 
into the conduction sea, we argue that the fixed point is ob- 
tained by fusing the tt/2 phase-shifted chain with an iso-spin- 
1 primary field. The spectrum obtained in this way agrees with 
the actual NRG results, see Table |l] [12] The approach to the 
fixed point is controlled by the leading irrelevant operator of 
dimension 7/5, which is the first descendant of the iso-spin-1 
primary field|ll] and represents the residual, actually attrac- 
tive, interaction among the conduction electrons once the im- 
purity has been absorbed. This operator leads to a singular 
local iso-spin susceptibility, xi ^ T^^/^, whose components 
are the compressibility and the pairing susceptibility in the 
S — L — Cooper channel. The fixed point has a residual en- 
tropy S{0) = 1/2 In [(75 + l)/(\/5 - 1)] and an impurity 
DOS p{0) = po/{l + V5), with p{e)~p{0) - |e|2/5.Jll This 
implies that the tunneling conductance approaches a fractional 
value at zero temperature, G(T = 0)/(2e'^/h) = 1/(1 + y/b), 
in a power-law fashion, G(T) - G(0) - T^/^. The fixed 
point is however unstable towards symmetry -breaking bound- 
ary terms which correspond in Table HI to physical operators 
with dimensions smaller than one. They include p-h and 
gauge symmetry breaking in the S = L = Copper-channel, 
i.e. the (/, S, L) — (1, 0, 0) operator, as well as quadrupolar 
distortions and spin-orbit coupling, the (/, S, L) = (0, 0, 2) 
and (0,1,1) operators, respectively. 

Since the two stable fixed points for J <C —Tk and 
J ^ Tk are essentially different, we expect an unstable 
fixed point in between, which we actually find by NRG for 
a J^, ^ —Tk. The NRG spectrum shows that the iso-spin 
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TABLE II: Same as TableUfor the unstable fixed point. This time 
however the states and operators are labelled by {Q, S, L, Z3) where 
Q refers to the U(l) charge sector and Z3 to the Potts sector. 
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symmetry is not restored at this point, unlike at the two sta- 
ble ones, which forces us to deal with independent charge 
U(l) and Potts Z3 sectors. The conformally invariant bound- 
ary conditions of a 3-state Potts model in two-dimensions are 
well characterized! 14| and include stable fixed boundary con- 
ditions, only one of the 3 states is allowed at the boundary, 
less stable mixed boundary conditions, two states are allowed, 
and lastly an unstable free boundary condition. We observe 
that the Kondo-screened phase and the non-Fermi liquid one 
can be identified respectively as the fixed and mixed bound- 
ary conditions in the Potts sector. Indeed the non-Fermi liquid 
fixed point spectrum can also be obtained fusing the Kondo- 
screened fixed point with the Z3 primary field e, in agree- 
ment with the boundary CFT of the 3-state Potts model ll4ll . 
The unstable fixed point should then corresponds to the free 
boundary condition. It has been shown| 14] that a proper de- 
scription of the free boundary condition in the 3-state Potts 
model requires a larger set of conformal towers which in- 
clude the so-called C-disorder fields llTsI of dimensions 1/40, 
1/8, 21/40 and 13/8. We actually find that the unstable fixed 
point is simply obtained by the Kondo-screened one upon 
fusion with the same 1/8 primary field that allows to turn 
the fixed into free boundary conditions in the 3-state Potts 
modell 14], see Table HII From the right panel in Table HII it 
turns out that the operator which moves away from the fixed 
point has quantum numbers (Q, L, Z3) = (0, 0, 0, Z) and 
dimension 2/3. This implies that the effective energy scale 
which controls the deviation from the unstable point behaves 
as r_ ^ |J — J*!'^. The most relevant symmetry -breaking op- 
erator has quantum numbers (Q, S, L, Z3) = (2, 0, 0, /) and 
corresponds to the S — L = Cooper channel. The corre- 
sponding pairing susceptibility diverges at zero temperature 
Xsc ^ T^^/^. The next relevant operators of dimension 2/3 
are the same quadrupolar distortions and spin-orbit coupling 
as for the stable non-Fermi-liquid fixed point. The residual en- 
tropy is S'(O) = 1/2 In 3 and the impurity DOS p(0) = po/2 
with /9(e) — p{Q) ^ — e^, implying a tunneling conductance 
G{Q)/{2e'/h) = 1/2. 




-J/Tk 

FIG. 1 : Qualitative phase diagram of the AIM as function of no and 
— J/Ta' for J < 0. Both screened and unscreened regions are Fermi 
liquids (FL). Screened/unscreened refers to the presence/absence of 
a Kondo resonance in the DOS. The critical line separating them is 
non Fermi liquid (NFL), as well as the line at no = 3 for —J/Tk 
above the critical value ~ 1. The arrow indicates the path along 
which the DOS in the right panel of Fig.|3|is calculated. 

Unlike the stable NFL fixed point, the unstable one survives 
p-h symmetry breaking. Therefore, although both regimes 
with ,/ ^ —Tk and J ^ Tk should have Fermi-liquid be- 
havior away from p-h symmetry, an unstable NFL fixed point 
is still expected to intrude between them. The NRG analy- 
sis confirms this CFT prediction not only slightly but also far 
away from p-h symmetry. For instance we still observe an un- 
stable fixed point for J* ~ —Tk when the average impurity 
occupancy no = 2. In this case for J ^ —Tk the impurity 
locks into a total singlet, S = L = Q, which decouples from 
the conduction sea. For J ^ Tk the impurity instead is in 
a state with S = L = I which, due to the (/, s) = (0, 0) 
term in Eq. Q, gets fully screened by four conduction elec- 
trons. These two phases turns out to be still separated by a 
critical point, which we find can be obtained by fusing the 
spectrum of the unscreened fixed point ( J ^ —Tk) with the 
dimension- 1/8 primary field of the extended Z3 sector The 
operator content at this unstable fixed point is exactly the same 
as in the case of ng — 3, thus showing that they are connected 
by a whole critical Une J* (no). We find that J* ^ —Tk 
for no G [2,3]. For no < 2, |J*| increases and merges 
into a mixed- valence critical point J* ~ — J7 for no ^ 1, 
as sketched qualitatively in Fig.^ 

Let us discuss now the behavior of the impurity DOS in 
the various phases, which is directly proportional to the SMT 
conductance. In Fig.|2]we draw the DOS across the unstable 
NFL fixed point for no = 3. We notice that the DOS displays 
two distinct features at low energy. There is a broad resonance 
which is smooth across the transition. On top of that, there is 
a narrow peak in the Kondo-screened phase which shrinks, 
disappears at the fixed point, and transforms in a narrow cusp- 
like pseudo-gap within the NFL phase. We argue| 16] that this 
dynamical behavior is controlled by two energy scales. One is 
the aforementioned scale which measures the deviation from 
the fixed point, r_ ^ | J — J* p, and sets the magnitude of the 
naiTow resonance as well as of the pseudo-gap, p{e) — p{{)) ^ 
p(0) |e/T_|2/5. The other is a high energy scale, T+ - | J*|, 
which sets the width of the broad resonance. 

Since it is unlikely that experimentally the system finds it- 
self near the unstable fixed point, in Fig.|3lwe draw the low- 
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FIG. 2: NRG impurity DOS at no = 3 with ?7 = 2 and 
po = 2 in units of half the bath bandwidth fl^ . The large fig- 
ure shows the DOS on a large scale at J = 0, where the Hub- 
bard side-bands and the Kondo resonance are both visible. In 
the inset we show the low-energy part across the unstable fixed 
point. The curves from the top one downwards correspond to J = 
0.0, -0.01, -0.017, -0.018, -0.02, -0.035, -0.043. 
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FIG. 3: Low-energy impurity DOS at no = 1, 2, 3. (a) J = 0. The 
peaks of the curves with decreasing no = 3, 2, 1 gradually move 
towards positive energies, (b) J — —0.075. no = 3 is the curve 
with the cusp at the chemical potential, no = 2 that with the large 
pseudo-gap and no = 1 the other one. 

energy part of the DOS with uq — 1,2,3 faraway from the 
unstable fixed point. In the figure we show the two cases: 3(a) 
J = 0, standard Kondo effect; 3(b) J < along a path like 
the one drawn in Fig. [2 Unlike the J = case (also rep- 
resentative of all positive J's), for J < the DOS turns out 
to be highly sensitive to the doping no, similarly to what is 
observed experimentally! 5|. 

We end by discussing the role of quadrupolar distortions 
which lower the icosahedral symmetry of isolated Ceo and 
presumably exist in real SMT devices. In their presence the 
unstable fixed point critical region is replaced by a crossover 
region, more or less smooth depending on the value of the dis- 
tortion. Yet the spectral properties inside the inverted Hund's 
rule region should still differ from a conventional perfectly- 
screened Kondo-like behavior, particularly when the sensitiv- 
ity to the doping is concerned. 

In conclusion we have shown that a realistic Anderson im- 
purity model of a Cqq molecule displays a phase diagram 
which includes non-Fermi-liquid phases and critical points 
with anomalous properties, most notably higher sensitivity to 
doping n and quadrupolar distortions than to magnetic field. 



especially for n > 1. These results indicate a possibly very 
interesting behavior of Ceo based SMT's. Provided Ceo could 
be endowed with more than one electron in a controllable 
way, e.g. adsorbed on proper metal substrates f?] or by al- 
kali doping|5|, fractional zero-bias anomalies with power- 
law temperature/voltage behavior might eventually show up. 
Moreover the enhancement of the pairing fluctuations which 
emerges out of our analysis suggests an anomalous behav- 
ior when superconducting leads are used, as for instance an 
increase of the zero-bias anomaly below the lead critical- 
temperature if much smaller than Tk- 

In view of recent speculations about the role of unsta- 
ble fixed points in single AIM's in the context of strongly- 
correlated models on large-coordination lattices fl?'], our re- 
sults may also have important implications for alkali-doped 
fullerenes. We will return to this in a later work. 
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